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Abstract 

In this research, a gamma-type, low-temperature differential (LTD) solar Stirling engine with two cylinders was modeled, constructed 
and primarily tested. A flat-plate solar collector was employed as an in-built heat source, thus the system design was based on a 
temperature difference of 80 °C. The principles of thermodynamics as well as Schmidt theory were adapted to use for modeling the 
engine. To simulate the system some computer programs were written to analyze the models and the optimized parameters of the engine 
design were determined. The optimized compression ratio was computed to be 12.5 for solar application according to the mean collector 
temperature of 100 °C and sink temperature of 20 °C. The corresponding theoretical efficiency of the engine for the mentioned designed 
parameters was calculated to be 0.012 for zero regenerator efficiency. Proposed engine dimensions are as follows: power piston stroke 
0.044 m, power piston diameter 0.13 m, displacer stroke 0.055 m and the displacer diameter 0.41 m. Finally, the engine was tested. The 
results indicated that at mean collector temperature of 110°C and sink temperature of 25 °C, the engine produced a maximum brake 
power of 0.27 W at 14rpm. The mean engine speed was about 30rpm at solar radiation intensity of 900W/m“ and without load. The 
indicated power was computed to be 1.2W at 30rpm. 

© 2007 Elsevier Ltd. All rights reserved. 

Keywords: Solar energy; Stirling engine; Flat-plate solar collector 


1. Introduction 

In view of energy management, solar energy is one of the 
most important renewable and non-depletable energy 
sources for generating power. There are several methods 
for converting solar heat into mechanical energy. One of 
these methods that theoretically associated with maximum 
efficiency is Stirling engine (or hot air engine). Stirling 
engine is a simple kind of external combustion engine. It is 
an old concept first proposed by Robert Stirling in 1816 
(UK, patent no. 4081). Engines based upon his invention 
were built in many forms and sizes. These engines offer the 
possibly of a high-efficiency engine with lower exhausted 
emissions in comparison with the internal combustion 
engine. Hot air engines are clean and efficient and ran 
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almost silently on any combustible material such as field 
waste and biomass. 

In 1970 and 1980s a huge amount of researches were 
conducted on Stirling engine for automobiles by companies 
such as General Motors and Ford . The main drawback is, 
the Stirling engine tends to run at a constant power setting 
that is not proper for automobiles. But this characteristic is 
perfect for applications such as water pumping. Studies 
about high-temperature Stirling engines have been exten¬ 
sively reported. 

In most models, the engines operate with a heater and 
cooler temperature about 923 and 338 K, respectively [1]. 

The thermal limit for the operation of high-temperature 
Stirling engines depends on the material used for its 
construction. Engine efficiency ranges from about 30-40% 
resulting in a typical temperature range of 923-1073 K, and 
normal operating speed range is from 2000 to 4000 rpm [2]. 
On the other hand, the low-temperature differential (LTD) 
Stirling engine is a type of Stirling engine that can operate 
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Nomenclature 

(P 

phase angle 



a 

displacer crank angle 

h 

convective heat transfer coefficient 

p 

power piston crank angle = (cp—ck) 

P 

instantaneous pressure 



m 

gas mass 

Subscripts and superscripts 

k 

constant 



V 

volume 

c 

cold outer space/plate 

T 

temperature 

C 

cold inner space 

Q 

heat 

0 

reference parameter 

r 

gas constant 

d 

displacer 

S 

exchanger plate area 

ds 

dead space 

W 

work 

H 

hot outer space/plate 

X 

working piston position 

h 

hot inner space 

y 

displacer position 

loss 

heat loss 

A 

difference 

P 

working piston 

7 

gas heat capacity ratio 

reg 

regenerator 

V 

efficiency 

* 

dimensionless 

V 

frequency 

** 

dimensionless 


dead space volume ratio = ( V xds j V d ) 

/ 

isothermal process 


compression ratio = ( V d V p ) = (SakoVOSp-^o) 

Pt 

point 

X 

internal temperature ratio = (T h /T c ) 




with relatively small temperature difference between hot 
and cold ends of the displacer cylinder [3]. The efficiency of 
a LTD Stirling engine is low in comparison with the high- 
temperature Stirling engines but it can be acceptable 
due to the availability of many low-temperature heat 
sources including solar energy which is inexpensive 
and safe. LTD engines may be of two designs. The first 
one uses single-crank operation where only the power 
piston is connected to the flywheel, called Ringbom engine. 
A short, large-diameter displacer rod in a precise-machined 
fitted guide has been used to replace the displacer 
connecting rod [3]. The other design is called a kinematic 
engine, where both the displacer and the power piston are 
connected to the flywheel through the cranks of the 
crankshaft [3]. Large amount of studies on LTD Stirling 
engines have been done during the recent century. Some of 
them are as follows: 

Kolin [4] tested a number of LTD Stirling engines, over a 
period of many years. In 1983, he presented a model that 
worked on a temperature difference between the hot and 
cold ends of the displacer cylinder as low as 15 °C. 

Senft [5] made an in-depth study of the Ringbom engine 
and its derivatives, including the LTD engine. His research 
in LTD Stirling engines resulted in a most interesting 
engine which had an ultra-low temperature difference of 
0.5 °C. It seems difficult to create any advancement better 
than this result. Senffis work in 1991 [6], showing the 
principle motivation for using Stirling and general heat 
engines, developed an engine operating with a temperature 
difference of 2 °C or lower. 

In 1993, Senft [4] described the design and testing of a 
small LTD Ringbom Stirling engine powered by a 60° 
conical reflector. He reported that the tested 60° conical 


reflector, producing hot end temperature of 93 °C under 
running conditions, worked very well. 

In 1997, Iwomoto et al. [7] compared the performance of 
a LTD Stirling engine with a high-temperature differential 
Stirling engine. Finally, they concluded that the LTD 
Stirling engine efficiency at its rated speed was approxi¬ 
mately 50% of the Carnot efficiency. 

In 2005, Kongtragool and Wongwises [8] investigated 
theoretically the power output of the gamma-configuration 
LTD Stirling engine. In this work, the power calculation 
was studied and discussed. 

In 2005, Kongtragool and Wongwises [9] presented the 
optimum absorber temperature of a once reflecting full- 
conical reflector for a LTD Stirling engine and a 
mathematical model for the overall efficiency of a solar- 
powered Stirling engine was developed. 

In 2006, Kongtragool and Wongwises [10] designed and 
constructed two single-acting, twin-power piston and four- 
power piston gamma-configuration LTD Stirling engines 
with heater temperature of about 589-771 K. The tested 
engines were equipped with regenerator to increase the 
thermal efficiency. 

Solar energy is a proper heat source for LTD Stirling 
engine. Solar collectors as a special kind of heat exchangers 
are usually employed to convert the solar radiation 
into thermal energy. If a hot air engine is equipped 
with a concentrated solar collector as a heat source, the 
working temperature would be high and heat efficiency 
would be effectively improved. But this system always 
has to trace the direction of sun in the sky. The tracer 
is a complex system and can waste part of the energy. 
Flat-plate solar collectors need no tracer system and 
can absorb diffused solar energy radiation as well. As 
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Fig. 1. A general type of LTD Stirling engine with gamma-configuration. 


described above there is no published paper on experi¬ 
mental results of the kinematic LTD solar Stirling engine 
powered by a flat-plate solar collector without regenerator. 
In this research, the simulation, construction and testing of 
a solar LTD Stirling engine with flat plate solar collector as 
an in-built heat source and without application of the 
regenerator were investigated. 


2. Mathematical model 


A schematic sketch of a general LTD Stirling engine with 
gamma configuration is shown in Fig. 1. The Stirling cycle 
is composed of two isothermal and two isochoric processes 
(Fig. 2). The procedures of finite dimension thermody¬ 
namics were employed to calculate the optimized compres¬ 
sion ratio according to the guidelines of Pierre Rochelle 
[11]. An optimization process, with work or efficiency as 
the objective function, is investigated. In the following 
description, the classical Schmidt assumption will be used 
except for the finite heat transfer between the working gas 
and the exchanger plates. 

The gas mass content inside the engine is expressed by 
the sum of gas mass contents in the expansion and 
compression spaces together with the gas mass inside the 
regenerator: 


m = m c + m h + m reg 

P 

xSp + (To — y)$d + JAds 


rTr 


+ 


P 


rT h 


+ fAdsI 


(l) 



where y 0 is the complete stroke of the displacer piston, P is 
the instantaneous pressure, m reg is the mass of the gas 
inside the regenerator, F hds and F cds are the dead space 
volumes on the hot and cold sides. Since this research has 
been devoted to design a LTD solar Stirling engine without 
regenerator, thus m reg was assumed to be zero in Eq. (1). 
Moreover, in the LTD Stirling engine, if a regenerator had 
been employed for improving the engine efficiency, the 
value of m reg would have been negligible compared to 
(m c +m h ). This can be easily concluded form Fig. 1 
according to the large volume of expansion and compres¬ 
sion spaces in comparison with the small unavoidable dead 
volume of the regenerator. 

To reduce the number of variables involved in the problem, 
dimensionless descriptions are defined as follows [11]: 


( 2 ) 
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* 

X = 


* 

y = 


X 

x 0 ’ 

(3) 

y 

V 

(4) 


where x 0 is the complete stroke of the power piston. The 
reference mass m 0 and the dimensionless mass nf are defined 
as follows: 


mg = 


PoV p 

rT r 


(5) 


nf = — = P*\[x* + (1 - y*K + Ccr t ] + -\?Z + £a h ] 
mo I t 


( 6 ) 


Assuming the same dead space volumes o — cr h = cr c and 
letting A = £(1 + 2d) and 5 = £((1 /t) — 1). Substituting ^4, 
5 and c in Eq. (6): 


m* = 5*[x* + ^ + B(y * + d)]. 

Other dimensionless descriptions are as follows: 


(7) 


P*(x*,/) = 


m 


❖ 


x* + A + 5(y* + o') 


( 8 ) 


mf(x*,y *) = '— = m 


m h 




m 0 


+ cr) 


t[x* + ^4 + 5(y* + o')] 


(9) 


* */n n * £0 + 


W hmax = W h(°> ] ) = m 


t[A + 5(1 + cr)] 


( 10 ) 


<min =ml{\,0) = m 








t[ 1 + ^4 + 5cr] 


(ii) 


The dimensionless mass transfer Am* between the volumes 
is defined as 


Am* = 


max min 


= m 


* 


(1 + ^ + cr) 


t [v4 + 5(1 + <r)][l + A + 5cr] 


( 12 ) 


To get the energy balance of the cycle, the dimensionless 
work W* and transferred heats during the isothermal 
processes are expressed as 


W* = 


w 


1 


Ai V 


o y p 


PoSpXo 


PS p dx = 


- (b P* dx* 


(13) 


02 = -fPj-J PdV c = i> P*( dx* -cd/). 


(14) 


where V c = xS p + (y 0 - + F cds . 


- = _E/ 

PoV p J 


0l = iTTr<f> PdV h = ^(f> P*dy* 


(15) 


where V h =yS d + Vm s - 


Solving Eqs. (13)—(15), we obtain 
W* = - £ P* dx* 


*pt2 


•pt4 

5*(x*, 0) dx* + / 5*(x*, l)dx* 

_^ptl ./pt3 


-0 


= —m 


* 


dx* 


\ x* A Bo 


+ 


dx* 


x* + A + 5(1 + cr) 


= —m* In(Z). 


( 16 ) 


where 


Z = Z(if, (T, T) = 


1 + 0/A + 5(cr + 1)) 


1 +(1 /A + Bo) 


Q h = -m* z- In (Z) 


B 


m* 

>'* C In (Z) 


(17) 


(18) 


O'* = 

t B 

To calculate the total heat absorbed by the working 
gas at the heat source, the heat lost in the regenerator is 
added to the transferred heat during isothermal process 
[ 11 ]: 


e reg = Am—y(r h - T c ) 

Am 1 

= PoV p -r(r- 1), 

mo y — 1 

£?loss Sreg(^ ^7reg)’ 

2loss a -.x ^ 


eLs = 


= Am*(r - 1)-^ = —m*Gk 

PoV P v 7 7-1 


reg? 


(19) 


( 20 ) 


( 21 ) 


where 


1 — rj 


reg 


^reg — 

y — 1 

and 


G = Gfe o', t) = 


5(1 + x4 + cr) 


M + B(l + d)][l + A + 5(7] 

The total heat absorbed by the working gas at the hot 
source and the total heat rejected to the cold sink are as 
follows: 


* 


e h 


h + eL 

* 


= m 


Ql + o 

—m 

(OA) -1) 

- In (Z) 


In (Z) — m Gk reg 


❖ 


L(i/r) — i Gktes 


( 22 ) 


* 


Cc 


= m 


o'l - efoss 

m* In (Z) 

WJ7) -1) 

In (Z) 


+ m Gk veg 


* 


_t((1/t) - 1) 


+ Gk 


reg 


(23) 
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Finally, the efficiency of the engine can be calculated as From Eq. (31): 


\W*\ _ —In(Z) 

~Qf~ [(t ln(Z)/(l - z)) + k reg G] 

T — 1 

X + kreg(l ~ T)(G/ In (Z)) 



The internal working gas temperature at expansion space 
is less than the temperature of the solar absorber and the 
air temperature inside the compression space is more than 
the cold sink. To calculate the real temperature of the 
working fluid inside the cylinder, heat transfer equations 
are used as follows [11]: 







Qc 

P 0 V p 


hS d T H T c / 
vPo Vp P H V 


where 



hS d Tu _ hT h 

vF’o-S'dJ'o vP oyo ’ 






Qt = 


Qx 


■h 


hSi (T h - T h ) = Nt(l- Th 


PoV p vP 0 V p . ’V T h, 

T c can be determined from Eqs. (23) and (26): 


T c = T c -m 


* 


In (Z) 


NS, Lt((1/t) - 1) 


+ Gk 


reg 


= Th 


T c 


— m 


* 


In (Z) 


-+ (--l \k 
X \x 


G 

reg hi(zj 


r H nb 

T h can be determined from Eqs. (22) and (28): 

Me 


= Tu-rri 


-ln(Z) z ^ 
l(i/t) _ i ^ regC " 


= Th < 1 — m 


* In (Z) 


NB 


— 1 — k 


1 


reg 


- 1 


ln(Z). 


(28) 


(29) 


(30) 


From Eqs. (29) and (30), the internal temperature ratio t 
can be obtained: 


- T *_ 1 — /[—l-((l/T)-l)g] 

To {TclT u )-l[{\lg + ((\/x)-\)Hy 

where H = H(k [eg , c, z, a) = k rc „(G/ In (Z)) and I = I(m*, 
N, cr, t) = m* (In ( Z)/NB ). 



Tc 

Th 


+ IH 


z[2I + 1 ]-IH 




If T c and T h assumed to be close together approximately, 
then t« 1 and it results that m* = 1 + A, Z = 1 
-(5/04(1 + A))), In (Z) = ~(B/(A(l+A))), G = ((5(1 + 
a+ <r))/0(l + d))), H = -A: reg (l + d + ff ), / s —(m*/N) 
(1 /(d(l + A))) and IH^m\k Tes /N)(((l +A + <r))/(A(l + A)). 
Thus H and / are independent of r [11]. The root of 
Eq. (32) is 


a +1) 


i + 71 + [miTc/Tu)+my (i + V 2 


(Tc/T h ) + IH 



Reintroducing the real value of x from Eq. (33) into the 
expression of W* and r\\ 


x — 1 

x + /c reg ( 1 - t)(G/ In (Z)) 

I — 1 ^ x — 1 

X + (1 - t )// _ x - (1 - x)k veg (\ + A + a) ’ 



IE* = In (Z) ^ -/;/ 




❖ 


B 


A(l+A) 


(i-d A)) 


m 


* 


(1+2<7)[1 + £(1+2<t)]’ 



The dimensionless work W* is reduced relative to the 
working piston swept volume but for the users a preferred 
reference volume is the overall volume, related to the 
engine bulk. This is equal to the working volume plus the 
displacer swept volume added to the dead space volumes. 
Then a new dimensionless work will be defined as [ii] 



| w\ 

5o(F P + V d + 2V ds ) 


W* 

1 + ^(1 + 2cr) 



2.1. Optimized compression ratio 

The first optimization goal concerning LTD Stirling 
engines was to determine the optimum compression ratio 
For given values of the parameters in Table 1, the 
compression ratio was incremented to find the optimum 
value of £ corresponding to maximum values of r\ and W**. 
This was done by a computer program to determine the 
optimized compression ratio (Fig. 3). 

In this research, an effort was devoted to design a solar 
Stirling engine powered by a flat plate solar collector 


Table 1 

Reference parameters 


r H ( K) 

T c ( K) 

A) (pa) 

cr 

v (revs ') 

h (Wm _2 K _1 ) 

To (mm) 

y 

Vreg 

373 

293 

100 000 

0.1 

0.5 

10 

55 

1.4 

0 
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Fig. 3. Efficiency and dimensionless secondary work versus compression ratio. 


without regenerator. The practical temperature of a 
flat-plat solar collector is around 100 °C or more. 
Therefore TH selected to be 100 °C and T c assumed to 
be about the wet bulb temperature of the environment. 
Since the regeneration process was ignored in this study, 
rj = 0. The convective heat transfer coefficient /z, 
assumed to be 10Wm“ K“ on both cold and hot 
plates [11]. 

It can be seen in Fig. 3, there are optimum values of the 
compression ratio. The optimized compression ratio 
assumed to be 12.5 to maximize the secondary work and 
efficiency. 

2.2. Internal temperature of expansion space and 
compression space 

Substituting the optimized compression ratio and para¬ 
meters of Table 1 and Eq. (33) into Eqs. (29) and (30), the 
internal temperature of expansion and compression spaces 
were computed. A computer program was written to solve 
the equations. The corresponding values of T h and T c were 
calculated as follows: 

T h = 337 K, T c = 328 K. 


2.3. Schmidt theory 

The Schmidt theory is one of the isothermal calculation 
methods for Stirling engines. This theory provides for 
harmonic motion of the reciprocating elements but retains 
the major assumptions of isothermal compression and 
expansion and of perfect regeneration. It, thus, retains 
highly idealized, but is certainly more realistic than the 
ideal Stirling cycle [1]. The assumption of simple-harmonic 
volume variation permits pressure, P , to be expressed as a 
function of crank angle, a, and leads to closed form 
solutions for work per cycle [12]. The Schmidt formula may 
be shown in various forms depending on the notations used 
and can be arranged for gamma-configuration Stirling 
engines [4]. Since the internal temperatures of expansion 
and compression spaces, considering the proposed model, 
were calculated in Section 2.2, the Schmidt theory with all 


its related assumptions can be employed to determine 
the theoretical output work and optimized phase angle 
of the engine. According to the mechanical configura¬ 
tion shown in Fig. 1, the Schmidt formula can be arranged 
as follows: 


V h = V hds + f (1 — cos (a)). 


(37) 


V c = V cds + f(l — cos (a — (p)) + f (1 + cos (a)), (38) 


m = — 


P 

r 


T c 
PV_ d 

2 rT, 


Vo v h 


n 


2a h 11 

-1-cos (a) + 2c< 

T T T 


+ \ — Tos(a — q>) + cos (a) + 1 

£ £ 


(39) 


2mrT c 
V d F(a) ’ 



where F( a) = (2<Th/x) + (1 jz )—(1 /t) cos(a) + 2cr c + (1 /£)— 

(I/C) cos (a— q>) + cos(a) + 1. 



/ 


P dV/j = mrT^ 



sin(a) 




W c = jPdVo 


mrT c f 2n sin(a — cp) 

m 


da 



mrT c f 2n sin(a — cp) 

m 


total = W h + W c = 


da. 


(43) 


fE t otai is the total work done per cycle. Using the results of 
Section 2.2 and the information of Table 1, the total work 
done per cycle and the optimized phase angle can be 
evaluated. 
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2.4. Total work done per cycle and optimized phase angle 


3. Engine design and construction 


A computer program was written to analyze the Eq. (43) 
and plot the total work done per cycle as a function of 
phase angle to calculate the optimum phase angle and total 
work (Fig. 4). 

In this research, the engine was designed with two 
cylinders, so the total work done per cycle is 
2x 1.18 = 2.36(J) theoretically without any regeneration 
process. 



Fig. 4. Total work done per cycle versus phase angle. 


Main engine design parameters are shown in Table 2 
according to mathematical model. The schematic diagram 
of the engine was shown in Fig. 5. It is designed with two 
separate cylinders and Gamma-configuration arrangement. 

The power cylinders were made of a steel pipe and the 
power pistons were made of an aluminum bar. The power 
pistons were turned to match the power cylinder bores. The 
clearance between the power piston and power cylinder 
was 0.01 mm x 15 mm thick Teflon pan was used to 
construct the power piston connecting rod. The solar 
absorber and the cold plate were made of 3 mm thick 
aluminum plates. The crank shaft was made of 12 mm thick 
steel shaft with eight Teflon cranks and crank pins. The 
crank shaft was supported by six self-align ball bearings 
that were located inside six wooden casing. The flywheel 
was attached to the middle part of the crankshaft. The 
flywheel was constructed form 15 mm thick Teflon pan 
with several 0.14 kg steel dead weights attached around it 
to decrease the speed fluctuations. 

4. Measurement apparatus 


Table 2 

Design parameters 
Displacer 

Bore x stroke (m) 0.41 x 0.055 

Power piston 

Bore x stroke (m) 0.13x0.044 

Phase angle 90° 

Flat-plate collector dimensions (m) 1 x 0.5 


The testing facilities are shown in Fig. 7. Since the engine 
speed was low, an accurate photo tachometer with 
+ 0.1 rpm accuracy was used to measure the engine speed. 
Cooling water was used in the cold sink. To measure the 
temperatures of the absorber and cold sink, some SMT-160 
thermal sensors were attached to the aluminum plates of 
heat exchangers. The accuracy of the temperature mea¬ 
surement was ±0.5 °C. To have a wide range of collector 
temperature, a flat reflector was attached to the wooden 


Flat-plate collector 


Displacer 

cylinder 

Displacer 

connecting 

rod 


Flywheel 
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a- 
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connecting 
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Frame 


Fig. 5. Solid model of the LTD solar Stirling engine without regenerator. 
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Fig. 6. View of the LTD solar Stirling engine without regenerator. 


casing of the collector to enhance the solar energy intensity 
on the absorber plate of the collector (Fig. 6). 

The solar intensity radiation was measured and recorded 
by a Casella solarimeter every 5 min (0-2000 Wm~ , 
+1 mmvWm - ). 

Keeping in mind that generally the rope brake dynam¬ 
ometer with a spring balance and loading weights is used to 
measure the engine torque at different engine speed and the 
engine torque can be determined from the difference of spring 
balance reading and loading weight. Since the speed of the 
LTD solar Stirling engine is clearly low compared to the speed 
of internal combustion engines, an attempt has been made to 
employ another dynamic method for evaluating the torque 
and brake power of this low-speed LTD Stirling engine. Some 
weights, a rope and a brake drum have been used to measure 
the engine torque and the brake power (Fig. 7). The weights 
were hanged around the rotary brake drum attached to 
crankshaft through the rope and then, rotating the brake 
drum would cause the hanging weights to be elevated. 
Therefore, the hanging weights would act as a braking force 
against the rotary brake drum and the engine torque can be 
computed by multiplying the downward force of the loading 
weights by the brake drum radius once the constant engine 
speed is displayed on the digital tachometer. The brake drum 
diameter was 0.02 m and it was directly attached to the axis of 
the engine crankshaft through a central hole inside it. 


Thermal sensor 



Fig. 7. Schematic diagram of the solar Stirling engine rig and testing 
facilities (r: brake drum radius, m : mass of the loading weights, T: engine 
torque, n: engine speed, g : gravitational acceleration). 


5. Results and discussion 


5.1. Indicated power 


The mean indicated power can be calculated approxi¬ 
mately from 

W total ^mean 


^Indicated — 


60 


(44) 


where W tota i is the total work done per cycle computed in 
Section 2.4. and n mean is the mean engine speed at mean 
collector temperature of 110°C and sink temperature of 
25 °C. n mean was measured to be about 30 rpm. The 
indicated power was calculated to be about 0.6 W for each 
cylinder and thus the total indicated power for two 
cylinders would be 1.2 W. 
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5.2. Engine performance 

The engine was constructed and primarily tested at 
Shiraz University from 1 to 6 August 2005. Fig. 8 shows 
the engine torque variations versus the engine speed. In this 
figure, the engine torque T is calculated by multiplying the 
brake drum radius r by the downward braking force mg as 
shown in Fig. 7. 

T = rmg, (45) 

where m is the mass of the loading weights and g is the 
gravitational acceleration. Fig. 9 shows the brake power 
versus engine speed. The brake power can be calculated 
from dynamics equations as [13] 

2nTn 

^Brake = ? (46) 

60 

where n is engine speed in rpm and T is the engine torque 
in Nm. 

0.22 
0.2 
0.18 
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Fig. 8. Engine torque versus engine speed. 
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Fig. 9. Brake power versus engine speed. 


According to Fig. 7, the load is gradually applied to the 
brake drum, and then the engine speed would be gradually 
reduced till a certain applied load would finally shut down 
the engine. The experiments were conducted during the 
time interval of 11.30-13:00o’clock, when the mean 
collector temperature was measured about 110°C with 
minimum fluctuations. The characteristics of torque and 
brake power variations represented in Figs. 8 and 9, are 
similar to the results of the research published by 
Kongtragool and Wongwises [10]. In these figures, it can 
be noted that the torque and the brake power decrease with 
an increase in engine speed. This reduction may be 
attributed to the lack of rather good heat transfer at 
higher engine speed as well as increasing in friction. 

5.3. Relationship of collector temperature and engine speed 
at no-load condition 

The no-load speed of the engine at various collector 
temperatures for two different sink temperatures is shown 
in Fig. 10. 

It can be illustrated that, the engine speed has been 
increased with an increase in the collector temperature. It is 
also observed that with decreasing the sink temperature, 
the engine speed is effectively increased. These improve¬ 
ments may be due to the high heat transfer potential at 
high-temperature gradients. 

5.4. The effect of regenerator application on the engine 
thermal efficiency 

In this research, the engine was designed without 
regenerator. An important question would be raised about 
the amount of efficiency reduction due to exclusion of the 
regenerator. In order to figure out the effect of regenerator 
application, a computer program was written to analyze 
Eq. (34). The effect of regenerator efficiency variations 



Fig. 10. No-load speed of the engine versus collector temperature at 
T c = 25 and 35 °C. 
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Fig. 11. Engine thermal efficiency versus compression ratio for different regenerator efficiencies. 


versus the theoretical thermal efficiency of the engine was 
investigated. All reference parameters assumed to be 
constant according to Table 1 except the regenerator 
efficiency. The results were shown in Fig. 11. 

It can be argued that the engine thermal efficiency was 
effectively improved for the regenerator efficiency of 1.0. 
Moreover, the engine thermal efficiency showed a dimin¬ 
ishing trend towards the regenerator efficiency of zero. It 
was also observed that the optimized compression ratio 
shifted to right due to an increase in the efficiency of 
regenerators. 

6. Conclusions 

In this research, the possibility of generating power 
from low temperature sources such as flat-plate solar 
collectors was investigated. A very simple two-cylinder 
LTD solar Stirling engine without regenerator was 
designed and primarily tested. Although the regenerator 
was omitted to determine the minimum possible output 
power and to simplify the engine structure, but results were 
incredible. 

The mean engine speed was measured to be 30rpm at 
collector temperature about 107-113 °C and sink tempera¬ 
ture of 25 °C which are approximately similar to the 
desirable reference parameters in Table 1. The procedure of 
finite dimension thermodynamics is suggested to predict 
the engine speed and the internal temperature of expansion 
and compression spaces. Therefore, the internal tempera¬ 
ture ratio can be calculated and used in Schmidt theory to 
compute a more reliable indicated power. The indicated 
power was calculated to be 1.2 W at 30 rpm and the brake 
power was measured to be about 0.1 W at 30 rpm. 

The application of an efficient regenerator is emphasized 
to increase the thermal efficiency. According to Fig. 11, for 


the regenerator efficiency of 1.0 at ideal conditions the 
engine thermal efficiency is 0.069 whereas for the regen¬ 
erator efficiency of zero, the engine thermal efficiency is 
about 0.0122. Therefore, by using an efficient regenerator, 
it is possible to increase the thermal efficiency six times 
more. 

Upon optimization, the optimal compression ratio was 
computed to be 12.5 with collector temperature of about 
100 °C, sink temperature of 20 °C and without any 
regeneration process. If an efficient regenerator was 
employed, the optimal compression ratio would be 
increased to 16 according to Fig. 11. 

The corresponding theoretical efficiency of the engine for 
the mentioned designed parameters was calculated to be 
0.012 for zero regenerator efficiency. 

Although the classical Schmidt theory is considered as 
one of the most widely used method in designing the 
Stirling engines, but its predictions seem to have limit¬ 
ations due to some unreal assumptions of this theory such 
as finite heat transfer between the working gas and the 
heat exchanger plates. In this study, an effort was 
conducted to improve this drawback in Schmidt theory 
by calculating the more accurate gas temperature inside 
the exchangers. This was done by employing appropriate 
heat transfer equations together with the proposed 
thermodynamics method to calculate more realistic work¬ 
ing gas temperatures. This finding results in a better 
estimation of indicated power of the LTD solar Stirling 
engine. 
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